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Abstract 

We systematically derive the Lax pair formulation for both discrete and continuum 
integrable classical theories with consistent boundary conditions. 
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1 The discrete case 



Quadratic Poisson structures first appeared as the well-known Sklyanin bracket pQ . A more 
general form, characterized by a pair of respectively skew symmetric and symmetric matrices 
(r, s) appeared in [2] in the formulation of consistent Poisson structures for non-ultralocal 
classical integrable field theories. Finally it was shown [3] that this was the natural quadratic 
form a la Sklyanin for a non-skew-symmetric r-matrix, reading: 



Ll, Lr 



r-r 77 , L Y L 2 



Lxir + r^Lt-Ltir + r^Lx. 



XI) 



A typical situation when one considers naturally a quadratic Poisson structure for the Lax 
matrix occurs when considering discrete or continuous integrable systems where the Lax 
matrix depends on either a discrete or a continuous variable; the Lax pair is thus associated 
to a point on the space- like lattice or continuous line [H [5] . Let us first examine the discrete 
case where one considers a finite set of Lax matrices L n labelled by n G N. 

Lax representation of classical dynamical evolution equations [6] is one key ingredient in 
the modern theory of classical integrable systems [7j-[TT] together with the associated notion 
of classical r-matrix [121 03] • Introduce the Lax pair (L, A) for discrete integrable models 
[9] (see also [H] for statistical systems), and the associated auxiliary problem (see e.g. [I]) 



t/v+i = L n if) n 
ip n = A n Vv- 



1.21 



From the latter equations one may immediately obtain the discrete zero curvature condition: 



L n — ln+i L n — L n A n . 



The monodromy matrix arises from the first equation (11.21) (see e.g. [1]) 

T a (X) = L aN (X)...L al (X) 



1.3) 



1-4) 



where index a denotes the auxiliary space, and the indices 1, . . . , N denote the sites of the 
one dimensional classical discrete model. 

Consider now a skew symmetric classical r-matrix which is a solution of the classical 
Yang-Baxter equation [121 13] 



n 2 (Ai - A 2 ), ri 3 (Ai) + r 23 (A 2 ) + ri 3 (Ai), r 23 (A 2 ) 
and let L satisfy the associated Sklyanin bracket 

L a (A), L b (n)\ = rab(X-fi), L a (X)L b (n) 



0, 



1.5) 



1.6) 
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It is then immediate that (11.41) also satisfies (11.61) . Use of the latter equation shows that the 
quantities trT(X) n provide charges in involution, that is 

hr T n (\), tr T m (^)} =0 (1.7) 

which again is trivial by virtue of (11.61) . In the simple SI2 case the only non trivial quantity 
is trT(X) = t(X), that is the usual "bulk" transfer matrix. In the bulk case in particular the 
zero curvature condition (II. 3p is realized by L n and [I]: 

A n (X, fi) = r\X) tr a {T a (N, n- X) r ab (X - fi) T a (n — 1, 1; A)} (1.8) 

where we define 

T a {n, m; A) = L an {X)L an ^ 1 {X) . . . L am (X), n > m. (1.9) 

We now generalize the procedure described in [I] for periodic boundary conditions to 
the case of generic integrable "boundary conditions". We propose a construction of two 
types of monodromy and transfer matrices, and associated Lax-type evolution equations, 
albeit incorporating a supplementary set of non-dynamical parameters encapsulated into a 
"reflection" matrix K(X). Any physical interpretation of the iT-matrix as a description of 
the "boundary properties" may not be appropriate in all cases. We should stress that this 
is the first time to our knowledge (see also [15]) that such an investigation is systematically 
undertaken. There are several related studies regarding particular examples of open spin 
chains [16] [17], however the derivation of the corresponding Lax pair is restricted to the 
Hamiltonian only and not to all associated integrals of motion. In this study we present a 
generic description independent of the choice of model, and we derive the Lax pair for each 
one of the entailed boundary integrals of motion. 

Subsequently we shall deal with two types of classical algebras, which are derived from 
two known types of consistent quantum boundary conditions. These boundary conditions 
are known as soliton preserving (SP), (see e.g. [H] [22]), and soliton non-preserving (SNP) 
[23] |24"[ [25] . SNP boundary conditions have been also introduced and studied for integrable 
quantum lattice systems [26]-[30]. From the algebraic perspective the two types of boundary 
conditions are associated with two distinct algebras, i.e. the reflection algebra [18] and the 
twisted Yangian respectively [31], [32] (see also relevant studies [251 12H1 [30], [33], [M]). It will 
be convenient for our purposes here to introduce some useful notation: 

r ab (X) = r ba (X) forSP, f ab (X)=rlf(X) for SNP 
r*M) = r ab (X) forSP, r* ab (X) = r* b a (-A) for SNP 
^(A) = r ba (X) forSP, r* ab (X) = r*j(-A) for SNP (1.10) 
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The two types of monodromy matrices will respectively represent the classical version of 
the reflection algebra M, and the twisted Yangian T written in the compact form: (see e.g. 
USE): 

{^(Ai), T 2 (A 2 )} = r 12 (Ai - A 2 )T 1 (A 1 )T 2 (A 2 ) - T 1 (A 1 )T 2 (A 2 )f 12 (A 1 - A 2 ) 
+T 1 (A 1 )f 1 * 2 (A 1 + A 2 )T 2 (A 2 ) - T 2 (A 2 )r 1 * 2 (A 1 + A 2 )T 1 (A 1 ) (1.11) 

where f, r*, f* are defined in (11.101) . In most cases, such as the A^_ ± r-matrices r 1 ^ 2 = r 21 
implying that in the SNP case r* ab = r* b . In the case of the Yangian r-matrix Tyi = ?"2i, 
hence all the expressions above may be written in a more symmetric form. 

In order to construct representations of (11. lip yielding a generating function of integrals of 
motion one now introduces c-number (non-dynamical) representations satisfying the purely 
algebraic condition (11 .lip since: 

{^(AO, tf±(A 2 )} = 0. (1.12) 

Taking now as T(A) any bulk monodromy matrix (11.41) built from local L matrices obeying 
(I1.6P and defining in addition 

f (A) = T _1 (-A) for SP, T (A) = T'(-A) for SNP. (1.13) 

one shows that representations of the corresponding algebras M, T, are given by the following 
expression see e.g. [T8l [35] : 

T(A) = T(A) K~(X) T(A). (1.14) 

For a detailed proof see e.g. [35] . 

Define now as generating function of the involutive quantities 

t(X) = tr{K + (\) T(A)}. (1.15) 

Due to ffTTTj) it is shown that [HI [35] 

|*(Ai), t(A 2 )}=0, Ax, A 2 GC. (1.16) 

Usually one considers the quantity In t(X) to get local integrals of motion, however for 
the examples we are going to examine here the expansion of t(A) is enough to provide the 
associated local quantities as will be transparent in the subsequent section. Finally one 
shows that time evolution of the local Lax matrix L n under generating Hamiltonian action 
of t(A) is given by: 

L n (/i) = A„ + i(A, fj) L n {ix) - A n (A, ji) L n (fi), (1.17) 
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where A n is the modified (boundary) quantity, 

A„(A, y) = tr a (^K^(X) T a (N,n; A) r ab (X - /i) T a (n - 1, 1; A) K~(X) f a (X) 

+ K+{X) T a (X) K-(X) f a (l,n-l; A) r* ab (X + //) f a (n, iV; A)) (1.18) 

where T(n, m; A) is defined in (jl.9p and 

T(m,n; A) = L am (A) . . .L an (A) n > m. (1-19) 

To prove (11.181) we need in addition to (11.61) one more fundamental relation i.e. 

{4(A), AO*)} = L 8 (A)^(A)L 6 (Ai) - L 6 (//)f: 6 (A + /*)L tt (A). (1.20) 

Taking into account (11.61) and the latter expressions we derive: 

{t(X), L b M} = tr a ^K+(X)T a (N,n + l;X)r ab (X-fi)T a (n,l;X) K-(X)f(X) 

+ K+(X) T a (X) K-{X)f a {l, n] X) f: b (A + /i) f a (n+l,N;X)yM 

- LM tr a (K+{X) T a (N,n; A) r ab (X - fi) T a (n - 1, 1; A) K~(X)f(X) 

+ K+(X) T (A) JC(A) f 8 (l,n-l;A) ^ 6 (A + /i) T a (n,iV;A)). (1.21) 

Expression (11.181) is readily extracted from (11.211) . 

Special care should be taken at the boundary points n = 1 and n = N + 1. Indeed we 
set: T(^,JV+1,A) = T(0, 1,A) = f(l,0,A) = T(N + 1, N, A) = 1 We should stress that the 
derivation of the boundary Lax pair is universal, namely the expressions (11.181) are generic 
and independent of the choice of L, r. Note that a different construction of representations 
of ( 11. lip was already given in a very general setting in [36]. It is related to the formulation 
of non-ultralocal integrable field theories on a lattice and extends the analysis of [2]. 

1.1 Example 

We shall now examine a simple example, i.e. the open generalized DST model, which may be 
seen as a lattice version of the generalized (vector) NLS model, (see also [37 1 1381 139| HU l [35] ) 
for further details). The open Toda chain will also be discussed as a limit of the DST model. 
We shall explicitly evaluate the "boundary" Lax pairs for the first integrals of motion. We 
focus here on the special case of the simplest rational non-dynamical r-matrices [JT] 

P N 
r(A) = - where P = V E {j <g> E j{ (1.22) 
A .^""i 
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P is the permutation operator, and (-Ejj)w = ^ik^jl- 

The Lax operator of the gl(j\f) DST model has the following form: 

JV-1 JV JV M 

L(A) = (A - x^X^En + + bj^ x^E Xj -J^X^E^ (1.23) 

3=1 3=2 3=2 3=2 

with Xn\ Xn being canonical variables. In [35] the first non-trivial integral of motion for 
the SNP case, choosing the simplest consistent value K ± = I was explicitly computed: 

N N JV-1 , JV-1 



H = - \ E N ' - 6 E E - i E ( x n^ + ft 8 *? 

n=l n=l j=l j=l 

JV-1 

where N n = ^xg ) X^" ) . (1.24) 

3=1 

Our aim is now to determine the modified Lax pair induced by the non-trivial integrable 
boundary conditions. We shall focus here on the case of SNP boundary conditions, basically 
because in the particular example we consider here such boundary conditions are technically 
easier to study. Taking into account (11.181) we explicitly derive the modified Lax pair for 
the generalized DST model with SNP boundary conditions. Indeed, after expanding (11.181) 
in powers of A -1 , and recalling (ll.22p we obtain the quantity associated to the Hamiltonian 

(HMD 

A( 2 ) = XE n - Xti^Eji + b ^~ 1] E lv n e {2, ... N} 

m 3+1 

A< 2) = A£ n - b^x^Eji + bJ2^t 1} Eij, 

3+1 3+1 

A# +1 = \E n - Ys^En + J2 X t 1] Eir (1-25) 

3+1 3+1 

It is worth stressing that in the sli case the SP and SNP boundary conditions coincide given 
that 

L _1 (-A) = V L*(-A) V, V = antid(l,...,l). (1.26) 

The equations of motion associated to the Hamiltonian (jl.24p may be readily extracted by 
virtue of 

L = \n {2 \ l\. (1.27) 

Alternatively the equations of motion may be derived from the zero curvature condition, 
which the modified Lax pair satisfies. It is clear that to each one of the higher local charges 
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a different quantify An is associated. Both equations (ll.27p . (II. 3p lead naturally to the same 
equations of motion, which for example in the sl 2 case read as: 

x n = x 2 n X n + bx n+1 , X n = -x n X 2 n - bX n _i, n £ {2, . . . N - 1} 
X\ = x\X\ + bx 2 , Xi = —X\X\ — bx\ 

xn = x 2 n Xn + Xn, Xn = —xnX^ — bXN-i- (1.28) 

Note that the Toda model [12] may be seen as an appropriate limit of the sl 2 DST model (see 
also [12]), and the corresponding boundary Hamiltonian, Lax pair and equations of motion 
are easily obtained from the expressions above (for more details see [35J). 



2 The continuous case 

Let us now recall the basic notions regarding the Lax pair and the zero curvature condition 
for a continuous integrable model following essentially [1]. Define \l/ as being a solution of 
the following set of equations (see e.g. [1]) 

— = U0M,A)* (2.1) 

— =V(x,t,\)V (2.2) 

U, V being in general n x n matrices with entries defined as functions of complex valued 
dynamical fields, their derivatives, and the spectral parameter A. The monodromy matrix 
from ( 12. ID may be written as: 

T(x,y,X) = Pexp\ [ U(x',t, X)dx'\. (2.3) 



The fact that T also satisfies equation (12. ID will be extensively used to get the relevant 
integrals of motion. Compatibility conditions of the two differential equations (12. ip . (12.21) 
lead to the zero curvature condition [8] [TO] 



U-V + 



U, V 



0, (2.4) 



giving rise to the corresponding classical equations of motion of the system under consider- 
ation. 

Hamiltonian formulation of the equations of motion is available again under the r-matrix 
approach. In this picture the underlying classical algebra is manifestly analogous to the 
quantum case. The existence of the Poisson structure for U realized by the classical r- 
matrix, satisfying the classical Yang-Baxter equation (11.51) . guarantees the integrability of 
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the classical system. Indeed assuming that the operator U satisfies the following ultralocal 
form of Poisson brackets 



U a (x, A), V b (y, ji) \ = r ab (\- fx), U a (x,\) +U b (y,n) 8(x-y), 



(2.5) 



then T(x,y, A) satisfies (jl.6p . and consequently one may readily show for a system on the 
full line: 



{lnt7-{r(x,i/,Ai)}, hxtr{T(x,y,\ 2 )}} = 



(2.6) 



i.e. the system is integrable, and the charges in involution -local integrals of motion- are 
obtained by expansion of the generating function hitr{T(x,y, A)}, based essentially on the 
fact that T satisfies (12. ip . 

Our aim here is to consider integrable models on the interval with consistent "boundary 
conditions", and derive rigorously the Lax pairs associated to the entailed boundary local 
integrals of motion as a continuous extension of the discrete case described previously. We 
briefly describe this process below for any classical integrable system on the interval. In 
this case one constructs a modified transition matrix T, based on Sklyanin's formulation 
and satisfying again the Poisson bracket algebras K. or T. To construct the generating 
function of the integrals of motion one also needs c-number representations of the algebra M. 
or T satisfying (II. lip . (11.121) . The modified transition matrices, realizing the corresponding 
algebras R, T are given by (I1.14p . where now T defined in (12.31) and T in fll . 13j) . The 
generating function of the involutive quantities is defined in (11. 15ft and one shows in this 
case as well: 



jt(x, y, t, Ai), t(x, y, t, A 2 ) i =0, Ai, A 2 G C. 



(2.7) 



In the case of open boundary conditions, exactly as in the discrete integrable models, we 
prove (for more details on the proof see |15j ) 



{2.1 



T a (0,-L,X), V b (x,fj) 
where we define 



' a (x,\,(j,)+ M a (x, A,//), U b (x, n) 



(x,A,/x) = T(0, x, X)r ab (\ - fx)T(x, -L, X)K~(X)f(0, -L, A) 
+ T(0, —L, X)K-(X)f(x, —L, X)r* ab (X + fi)f (0, x, A). 



Finally bearing in mind the definition of t(X) ( 11.151) . and (12. 8ft we conclude with: 

dY(x, A, y) 



In t(X), U(x,n) 



dx 



+ 



V(x,A,/i), U(x,fi) 



(2.9) 



(2.10) 
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where 

V(x,\,fi)=r 1 (\) tr a (K + {X) M a (x, A, . (2.11) 

As in the discrete case particular attention should be paid to the boundary points x = 
0, —L. Indeed, for these two points one has to simply take into account that T(x, x, A) = 
T(x, x, A) = I. Moreover, the expressions derived in (12.91) . 02. 1 lj) are universal, that is 
independent of the choice of model. 



2.1 Example 

We shall now examine a particular example associated to the rational r- matrix 01 .221) . that 
is the glj\f NLS model. Although in [35] an extensive analysis for both types of boundary 
conditions is presented, here we shall focus on the simplest diagonal (K ± = I) boundary 
conditions. The Lax pair is given by the following expressions [H EH] : 



U = U + AUi, V = V + AVi + A 2 V 2 (2.12) 



where 



j Af-1 Af-1 



t=l i=l 

JV-l Af-1 



i, j=l i=l 

Vi = -U , V 2 = -Ux (2.13) 



and ipi, tpj satisfjQ: 

|^(x), 4)j(y)\ = lA(x), $j(y)\ = 0, |^(x), $j(y)\ = 5 i3 - K x ~ v)- ( 2 - 15 ) 
The boundary Hamiltonian for the generalized NLS model may be expressed as (see 

rO Af-1 N-l 



n = I dx ^(k|^(x)| 2 ^|^-(x)| 2 + ^i 

J ~ L i=l j=l 



Af-1 Af-1 

E ($(0)&(0) + V*(0)#(0)) + Y, {^{-L)U-^+i^-m{-L)){2.lQ) 



i=l i=l 



4 The Poisson structure for the generalized NLS model is defined as: 

" L / SA SB SA SB 



\A, B}-i^J^dx{-^-^ _____ __j (2.14) 



s 



One sees here that the K- matrix indeed contributes as a genuine boundary effect. The 
Hamiltonian, obtained as one of the charges in involution (see e.g. [35j for further details) 
provides the classical equations of motion by virtue of: 

*M*,t) _ r w(0i _ L)j MXtt) \ = r w(0| _ L)| 



dt i n rtv ' y r at 

-L < x < 0. (2.17) 

Indeed considering the Hamiltonian H, we end up with the following set of equations with 
Dirichlet type boundary conditions 

^(0) = ^(-L) = z G {1, . . . , jV — 1}. (2.18) 

For a detailed and quite exhaustive analysis of the various integrable boundary conditions 
of the NLS model see |35|. 



As mentioned our ultimate goal here is to derive the boundary Lax pair, in particular the 
V operator. Hereafter we shall focus on the SP case with the simplest boundary conditions 
i.e. K ± = I. For any glj^ r- matrix we may expand (12. lip , taking also into account (ll.22p . in 
powers of A -1 (we refer the interested reader to [35l [T5] for technical details) and we obtain 
N^\x, A) -the bulk part- coincides with V defined in 02.12p . (12.131) . and for the boundary 
points Xb 6 {0, — L} in particular: 

^ 2 Af- 1 Af-1 

Y {3) (x b , A) = -— ( E n ~ + i Yl $i{xb)il)j{xb)Eij - i (fi(x b )EiH- - ^(x^E^ . 



i=l ij=l i , j = 1 



(2.19) 



We may alternatively rewrite the latter formula as: 

V (3) (x b , A) = Y{x b , A) \^i{x b )\ 2 E NN + \J2$i(xb)E w + ^{x b )E Mi ). (2.20) 

i=l i=l 

The last two terms additional to V (12.121) . (I2.13P are due to the non-trivial boundary condi- 
tions; of course more complicated boundary conditions would lead to more intricate modifi- 
cations of the Lax operator V. 

Acknowledgments: We wish to thank the organizers of RAQIS'07 where this work was 
presented. J. A. thanks INFN Bologna, and University of Bologna for hospitality. This work 
was supported by INFN, Bologna section, through grant TO 12. 
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